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Attractive Casimir Forces in a Closed Geometry 
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We study the Casimir force acting on a conducting piston with arbitrary cross section. We find the 
exact solution for a rectangular cross section and the first three terms in the asymptotic expansion 
for small height to width ratio when the cross section is arbitrary. Though weakened by the presence 
of the walls, the Casimir force turns out to be always attractive. Claims of repulsive Casimir forces 
for related configurations, like the cube, are invalidated by cutoff dependence. 
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In 1948 Casimir predicted a force between conducting 
plates due to quantum electromagnetic fluctuations [lj. 
Casimir forces have now been measured precisely 0, and 
agree with theory to a few percent accuracy. As minia- 
turization continues, static friction due to the Casimir 
attraction between components in microelectromechani- 
cal systems (MEMS) has become a problem of increasing 
concern^. Since the Casimir force is strongly geome- 
try dependent 0, the question arises whether it can be 
made repulsive by arranging conductors appropriately. It 
has been claimed that the parallelepiped of sides a, 6, b 
has positive Casimir energy for aspect ratio in the range 
0.408 < a/b < 3.48 HQ. If so, a repulsive force would 
occur for a/b > 0.785, as a is varied at fixed b. 

However this example, and indeed all claims of repul- 
sive Casimir forces, require elastic deformations of single 
bodies (a rectangle or a circle in 2-dimensions, a paral- 
lelepiped, a cylinder or a sphere in 3-dimensions, typi- 
cally) treated as perfect, often infinitely thin conductors. 
It has been shown that, when the conductors are modeled 
more realistically, the Casimir energies associated with 
deformation actually depend strongly on material prop- 
erties such as the plasma frequency and skin depth|7J, 
and diverge in the perfect metal limit where these cutoffs 
are ignored. In contrast, the forces between rigid bodies 
remain finite in that limit. Recently Cavalcanti 8] intro- 
duced a modification of the rectangle — the "Casimir pis- 
ton" — that is demonstrably free of cutoff dependence 0. 
The 2-dimensional piston studied in Ref. |8| consists of a 
single rectangle divided in two by a partition (the "pis- 
ton" ) . Cavalcanti calculated the Casimir force on the pis- 
ton due to fluctuations of a scalar field obeying Dirichlet 
boundary conditions on all surfaces. He found that the 
force on the piston is always attractive, although sub- 
stantially weakened with respect to parallel lines. 

In this Letter and the coming paper [TT) we consider 
the 3-dimensional piston, for both scalar and electromag- 
netic (EM) fields. We keep careful track of possible cut- 
off dependences and show that they cancel for the piston 
configuration0. First we give the exact result for pis- 



tons of rectangular cross section, expanding the result in 
powers of the piston-base separation, a, and identify the 
terms with specific optical paths Next we consider 
pistons of arbitrary cross section, where an expansion for 
small a can be derived. Finally we show that the force 
on the piston is always attractive for a rectangular cross 
section and argue that this is true for any cross section. 

The Casimir energy for a domain V can be defined by 
(H = c = 1), E(A) = |£™u; m (X>) - E Q , where u m (V) 
are eigenfrequencies in the domain T>, E$ is the energy 
of the vacuum, and A is a cutoff. The Casimir energy 
can be related to the Green's function, G(x',x, £), (here 
£ = lo 2 ) that obeys (-V' 2 -£)G(x', x,£) = <5 3 (x'-x) and 
the appropriate boundary condition, through: E(A) = 

\l K2 d£^f£p{£) where p{£) cx 3 Jdx G(x, x, £). We 
find it useful to expand G(x',x, £) in a series of "op- 
tical paths" as in Ref. 01 • The sum over optical paths 
isolates cutoff dependences, highlights the most impor- 
tant contributions, and is exact for rectilinear geometries, 
provided paths touching edges and corners are properly 
included [llj]. (We explicitly checked the results for the 
parallelepiped by other methods.) The Green's function 
is given by a discrete sum over straight line paths, p r 
(some of which are shown in Fig. from x to x' (of 
length l Pr (x' , x) ) , undergoing specular reflection r times 
from surfaces (eg. Fig.^a), (c), and (f)), edges (Fig.^b) 
and (g)), and corners (Fig. ^h)). Each path is weighted 
by a phase, 4>(p r ) = r) ns+nc , where n s and n c are the 
number of surface and corner reflections on the path and 
rj = — 1(+1) for Dirichlet (Neumann) boundary condi- 
tions (BC) . The cutoff independent, a-dependent part of 
the Casimir energy in a rectilinear geometry is 



(1) 



here i Pl .(x) is the coincidence limit (x' — > x) of path 
lengths. The cutoff dependent contributions (in the sense 
of |9j) are isolated in the one-reflection term, and has 
been omitted from eq. Some one- reflect ion paths are 
shown in Fig. ^a) and^Jb). 
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From considering the cutoff dependent contributions of 
one reflection paths in Region I (see Fig.^J, we obtain (lT) 
the following Casimir energy as A 
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where S — 2(ab + be + ac) is the surface area and 
L = 4(a + b + c) is the perimeter length. The cut- 
off dependent terms agree exactly with the predictions 
of Balian and Bloch's theory of the density of states in 
bounded domains [13j. For the geometries considered in 
Refs. [H, these divergences are not cancelled by other 
contributions eg. from the exterior (discussed below for 
the electromagnetic case), but are instead dropped as 
part of a physically unmotivated "renormalization" pro- 
cess. In reality, deformations of a single parallelepiped 
that change S or L cannot be calculated independently 
of the material properties of the metal that determine the 
cutoff. However, when the contribution of the upper part 
of the piston (Region II) is included, the a-dependence 
disappears from the cutoff dependent terms, leaving a 
cutoff independent force in the limit A — > oo. For a scalar 
field obeying Dirichlet (-q = —1) or Neumann (77 = +1) 
BC, and setting b = c for convenience, the force is 
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where A = b 2 and P = Ab are the area and perimeter of 
the base, respectively. We have expanded in a <C b = c 
and taken the total height of the box h — > 00, the con- 
figuration most relevant to potential experiments. Com- 
plete results will be presented in Ref. [ll]]. Here J n = 
Z 2 (l, 1; 4) + 27T77£(3), where Z m (xi,X2, ■ ■ ■ , x m ; n) is the 
Epstein Zeta functional and C,(n) is the Riemann Zeta 
function, such that J_i as —1.5259 and J+i s» 13.579. 

The terms in eq. (0 can be understood through the 
optical paths associated with them. The most important 
contributions have the shortest path length, l Pr (x) . The 
leading contribution, ~ 1/a 4 , comes from paths that re- 
flect only off both the base and the piston, beginning with 
the 2 reflection path, Fig. ^c), and continuing through 
4 (Fig. IHd)), 6, 8, • • • reflections. They sum up to the 
well-known parallel plate result, the first term in eq. 
The first correction, ~ 1/a 3 , comes from closed but not 
periodic paths that additionally reflect once from a wall 
of the piston, beginning with the three reflection path, 
Fig. ID^e), and continuing through 5 (Fig. Hff)), 7, etc. 
Their sum is related to the Casimir force between lines in 
two dimensions |l if. The next correction, ~ 1/a 2 , comes 
from paths that reflect off edges, beginning with the three 
reflection path, Fig. Cfg). It is related to the Casimir 
force between points on a line. The a-independent fourth 
term in eq. iJSJ arises from paths in Region II that do not 




FIG. 1: Optical paths in a 3-dimensional box of size hxbx c 
with a partition at height a. The origin and end points (which 
actually coincide) of some paths are separated so the path can 
be visualized. See the text for further description. 



reflect from the piston (see Fig.P^i)). Note that the con- 
stant term gives a repulsive (attractive) contribution for 
Dirichlet (Neumann) BC. The remaining terms are expo- 
nentially small as a — » due to an effective mass for the 
corresponding paths, which we explain in detail later. 

We now turn to the case of electromagnetism, which is 
of most physical interest. To construct the Casimir en- 
ergy, we must determine the frequency spectrum. Inside 
the cavity E and B obey the wave equation, V 2 (E, B) — 
<9 2 (E,B) = 0, supplemented by conducting boundary 
conditions at the walls, n x E = and n • B = 0. In 
Region I, the solution for E is, 



E r 



A x cos(n7ra/a) sin(m7ry/6) sin(^7rz/c)e 
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Here A x<y>z are constants and 



and similarly for E y , E 
n, m, / are integers. The magnetic field is obtained from 
B = -iV x E/w. The transversality of E, V • E = 0, 
restricts the degeneracy of the modes: If all integers are 
non-zero, there is a twofold degeneracy. If one of the 
integers is zero, there is a single solution. If two or more 
of the integers are zero, the mode is forbidden. Taking 
note of this, we may express the electromagnetic Casimir 
energy (Eq) in Region I in terms of the Dirichlet (Ed) 
and Neumann (En) energies, 
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Eh (A) = E D (A) + E N (A) -J^Ei (* , A) , (5) 

»=i 

where E\ (d% , A) is the Casimir energy in 1-dimension 
with d\ = a, d,2 — b, (I3 = c. The cutoff, A, reminds us 
that these Casimir energies are all cutoff dependent. We 
have already constructed the first two terms in eq. 
and E 1 (d i ,A) = d i A 2 /27r - ((2)/(47rd l ) is well known. 
Note that any terms that involve r/ in E^ cancel upon 
summation. This includes the surface divergence ~ SA 3 , 
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but the perimeter divergence ~ LA 2 remains. In the liter- 
ature, the possibility that this term is cancelled by a sim- 
ilar contribution from the exterior of the parallelepiped 
has been suggested. If the conducting walls have zero 
thickness and yet remain a perfect conductor, this can- 
cellation may occur. However, we believe this to be un- 
physical and contrived. Furthermore, accounting for the 
energy from the exterior of the parallelepiped is essential 
to determine its geometry dependence. In the piston ge- 
ometry, all cutoff dependence drops out, and exterior 
energies are handled unambiguously. 

Using the methods outlined for the scalar field together 
with the decomposition in eq. (jSJ) , we obtain an exact 
expression for the electromagnetic Casimir force on the 
piston with 6 = c in the limit ft — > oo, 
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with fmn(x) = 7rV m 2 +n 2 x and {m, n} S 1? \ {0,0}. 
The sum in the last term is highly convergent and ex- 
ponentially small as a — > 0. Temperature corrections, 
calculated using the techniques in Ref. 0], wm be pre- 
sented in Ref. [ll|. The first term in eq. © is Casimir's 
famous result for the force between two parallel plates 
of area A separated by a distance a Note that the 
perimeter term ~ P/o? in eq. © cancelled between the 
two electromagnetic polarizations. The remaining first 
correction, an edge effect, ~ 1/a 2 , provides a repulsive 
contribution. The coefficient of the a-independent part 
is J c = J-x + J+i = 2Z 2 (1, 1; 4) w 12.053. 

Before discussing the implications of this result, we ex- 
tend it to pistons of arbitrary cross section. Solutions to 
the wave equation in a piston of arbitrary cross section 
factor into the product of solutions on the line [0,a] and 
solutions in 5cl 2 , the cross section of the piston. The 
density of states is therefore a convolution of the two 
dimensional density, p${£), and the density on a line, 
Pi{£). Pi(£) is trivial and the smoothed ps{£) can be 
found in Ref. 0] for Dirichlet and Neumann BC. A de- 
composition of the EM case into Dirichlet and Neumann 
sub- pro blems, similar to eq. (JSJ, exists (for details, see 
Ref. ) - The results for scalar and EM fields are 



F n = 



F c = - 



3C(4) + 



167r 2 a 4 ' 327ra 3 Aita 2 
3C(4)_ A , C(2)(1-2 X ) 



0(1). 



8n 2 a 4 Ana 2 
Here \ depends on the 2-dimensional cross section, as 
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where on is the interior angle of each sharp corner and 
/c(7j) is the curvature of each smooth section described 



by the curve jj. For example, x — 1/4 f° r a rectangle 
and % = 1/6 for all smooth shapes. Substituting \ = 1/4 
we find the first three terms of eq. and the first two 
terms of eq. (jSJ). The a-independent terms in eqs. © 
and ©, denoted by — J/ A, are difficult to recover for 
arbitrary cross sections with the methods used here. 

We now discuss the implications of our results for issues 
of attraction versus repulsion. The term proportional to 
— Jc/A in eq. © plays an important part. To under- 
stand its origin consider the case ft » i » a. In this 
regime the cutoff independent contribution to Fc from 
Region I contains no a-independent term. The terms that 
scale like 1/a 4 and 1/a 2 remain. They come from paths 
whose lengths vanish as a — > 0. The exponentially small 
terms come from paths that hit both sides of the cavity. 
Their contributions vanish like e ~ 27rb / a , since the mini- 
mum wave number in the vertical direction, ~ 1/a, acts 
like a mass for horizontal propagation. The —Jc/A term 
in eq. I© is instead the sole cutoff independent relic of 
Region II when ft > fc> a. It comes entirely from opti- 
cal paths lying in the horizontal plane, such as Fig.^i). 
These contributions are extensive in the vertical direc- 
tion, and thus depend on a as ~ —(ft — a)/b 2 . This re- 
pulsion between the top and bottom plates of Region II 
leads to the a-independent attraction in eq. © . All other 
optical paths in Region II give contributions that vanish 
as ft — > oo. Now consider Region I when a> ft. Following 
the same logic, an extensive contribution to the energy, 
~ —a/6 2 , now arises from horizontal paths in Region I. 
This cancels the contribution from horizontal paths in 
Region II leaving a Casimir force, Fc, that approaches 
zero (from below) exponentially fast as a/b — > oo. This 
can easily be understood: In this limit the horizontal 
wave number, ~ 1/6, now plays the role of a mass for 
propagation in the vertical direction, which resembles a 
1-dimensional system. This cancellation occurs for any 
cross sectional shape, differing only by the value of Jc- 

Had we not formed the piston, but only considered the 
finite part of the energy in the "box" defined by Region 
I alone, the term — Jc/(32tt 2 A) would have been absent 
from the expansion equivalent to eq. ©■ As a result 
Fbox — > + Jc/(327r 2 /l) for a ^> 6. In particular, one can 
show Fbox > 0, for a/6 > 0.785. This explains the claim 
of repulsion that has appeared in the literature: Without 
Region II (or some open region that allows rigid motion 
of the partition) the Casimir energy of the parallelepiped 
is, in fact, cutoff dependent. If the cutoff dependence is 
somehow ignored, a repulsive force at large a/6 remains 
as an artifact. In Fig. [3 the force is plotted normalized 
to the parallel plates result F\\ = -3C(4)A/(8?r 2 a 4 ). The 
solid line gives our exact result, and the dashed lines 
are the results from successively including each of the 
finite corrections of Eq. |JBJ • There is good agreement for 
a < 0.36 with the first correction, and for a < 0.76 upon 
inclusion of the second (constant) term. 

In summary, we have computed the Casimir force act- 
ing on a rectangular piston (Fig.QJ exactly for both scalar 
and EM fields — one of a rather small set of exact results 



4 




0.2 0.4 0.6 0.8 1 .,1.2 1.4 

a/b 



FIG. 2: The force on a 3-dimensional piston Fc with a square 
cross section due to EM quantum fluctuations, as a function 
of a/b (see Fig. 0, normalized to the parallel plates result 
). The exact expression is the thick 
solid line, and the first and second corrections are the dashed 
lines. The thin solid line is the force -Fb ox that comes from 
differentiating the finite part of the energy in a single box. 



for physically interesting geometries. We make quantita- 
tive predictions for the forces in this physically realiz- 
able system (Fig. [5J , which are found to be attractive for 
any rectangular cross section. Casimir force calculations 
for isolated parallelepipeds that appear to give repulsion 
have ignored cutoff dependence The simplest way to 



cancel the cutoff dependence is to introduce Region II, 
which is found to also cancel a cutoff independent part of 
the energy. This changes the repulsion into attraction, 
and is expected to occur for any cross section. 

The corrections to the leading parallel plate force can 
be understood by analyzing an optical path representa- 
tion of the Green's function. For a scalar field the first 
correction is one power of a suppressed, and comes from 
paths that reflect only once from the piston's sides. It is 
repulsive for Dirichlct and attractive for Neumann BC. 
For the electromagnetic field this correction cancels be- 
tween Dirichlet and Neumann modes. The first non-zero 
correction is two powers of a suppressed and comes from 
paths that reflect once off an edge. It is repulsive. 

We also obtain the leading correction to the parallel 
plate result for arbitrary shape of the cross section. This 
correction is quite sensitive to the geometry of the base; it 
is the same for all smooth curves, but distinct for shapes 
with sharp edges. Note that the expansion parameter a 
provides a length scale for smoothing irregularities of the 
shape at small scales. 
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